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SUMMARY 
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Vi ' \ 

The equilibrium temperature distribution, thermal emission 
spectrum, and limb function of a nongrey N 2 -CO 2 -H 2 O model atmos- i 
phere have been computed for a large variety of atmospheric para- \ 
meters which are appropriate for the atmosphere of Venus in the \ 
vicinity of the cloud tops. Near infrared absorption by CO 2 of 
direct and diffusely reflected solar radiation as well as nonblack 
cloud emissivities have been included. 

The important effects of the near infrared solar heating by 
C0 2 absorption, the water vapor amount, and the cloud top pressure 
level on the computed results are demonstrated. An application of 1 

this study to several observations of Venus yield the conclusion \ 

\ 

that atmospheric absorption by CO 2 cannot account for the limb 
darkening observed in the 8-13H region. Other conclusions 
resulting from a study of the water vapor abundance and thermal 
emission spectrum depend on the cloud composition. If the clouds 
are ice crystals, their elevation is at the 0.2-0. 4 atm pressure 
level. If their composition is nonaqueous, a pressure level of 
0.4-0. 5 atm may be assigned to the cloud tops. These results are 
based on the assumption of a single reflecting cloud layer, an 
assumed CO 2 concentration of 5% and a brightness temperature 
of 225°K for the cloud tops in the 8-13n spectral interval. 

Several observations are suggested in order to clarify the 
relative importance of the effects of various phenomena on the 
atmospheric and cloud structure of Venus. 

/ 
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I . INTRODUCTION 


In a previous study, the authors (Reference 1) considered the 
structure of a planetary atmosphere in a state of radiative equili- 
brium. The radiative balance of a nitrogen-carbon dioxide atmos- 
phere was determined to yield a number of temperature distributions. 
In that non-grey study, the absorption spectrum of CO 2 was divided 
into 18 intervals. Absorption coefficients were derived for each 
of these intervals from the calculations of Stull, Wyatt and Plass 
(Reference 2) and from the laboratory data of Burch, Gryvnak, 
Singleton, France and Williams (Reference 3). The strong line 
approximation (Reference 4, for example) was found to be valid for 
the range of pressures, temperatures, and concentrations under 
consideration. Furthermore, a plane parallel atmosphere consisting 
of 34 layers of constant mixing ratio, in local thermodynamic 
equilibrium was assumed. A lower surface with a specified tem- 
perature and emissivity and solar radiation incident from above 
completed the boundary conditions. Radiative balance was obtained 
by an iterative calculation on a digital computer. 

A large number of solutions were found for various CO2 concen- 
trations, solar zenith angles, and cloud top temperatures, 
emissivities and pressures. Each solution yielded a temperature 
profile, the associated spectral intensities, and the limb 
functions . 

The primary results of these studies, as applied to a portion 
of the Venus atmosphere, indicated that near-infrared absorption 
of solar radiation by CO 2 causes considerable atmospheric heating. 
Temperatures above the cloud layer may differ by as much as 70°K 
between the subsolar point and the night side. The cloud top 
pressure level also produced significant effects while the 
variation of the CO 2 concentration and the cloud top emissivity 
was less pronounced. 

A number of new features were introduced into these calcula- 
tions to improve and generalize their applicability. The changes 
were made desirable partly by the results obtained from that study 
and partly by new observational data. 

Some previously obtained solutions showed large temperature 
gradients near the cloud top surface. In some cases, the adiabatic 
lapse rate was exceeded and, according to the Schwarzschild 
stability criterion, convective transport of energy would take place. 
A more realistic treatment should therefore include convection, 
and hence the calculations were modified to incorporate these 
effects . 
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In applying these calculations to the atmosphere of Venus, a 
second modification was found necessary, as a result of recent 
observations which confirmed the presence of water vapor. Hence, 
an analysis of absorption laws, similar to that carried out for CO 2 , 
was applied to the water vapor spectrum and absorption coefficients 
were derived. 

The final change concerns the method of obtaining the tem- 
perature distribution. This change is minor although it is expected 
to yield somewhat more accurate temperatures. 

The above modifications permit a study of the following items: 

1) The effects of water vapor on the equilibrium 
temperature profiles, 

2) The effects of atmospheric parameters on the 
convective motions arising from convectively 
unstable solutions, 

3) The dependence of the thermal emission spectra 
and the limb functions on the atmospheric 
parameters , 

4) A synthesis of the above analysis into solutions 
which may be compared with planetary observations. 

Before these items are illustrated, a brief discussion of the 
various modifications applied to the calculations will be presented. 

II. NEW FEATURES OF THE CALCULATIONS 

A. Water Vapor Absorption 

The absorption spectrum of water vapor consists of three 
intervals of interest: the near infrared spectrum from 1 to 5|j 
with numerous bands, the 6.3p band extending from 5 to 8|i, and the 
pure rotation band covering the interval from 10 to 1000|i. 

For this analysis, the 6.3|a band is divided into two and 
the rotation band into 4 intervals as shown in Table I. Although 
the rotation band extends below 20|j and overlaps the CO 2 absorption, 
it is ignored in this range with only the CO 2 absorption taken 
into account. This simplification appears justified, since the 
optical thickness of H 2 O in the 10-20^i region is at least one 
order of magnitude smaller than that of CO 2 for the concentrations 
expected on Venus. For the same reason, the near infrared 
absorption of H 2 O was ignored. 

Water vapor absorption data are required between 150-300OK 
and for the pressure range 1 to 10“ 3 atm. The absorption spectrum 
of water vapor has been studied in detail by a number of investi- 
gators. A convenient summary has been given by Goody (Reference 5). 
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TABLE I 


! 

Spectral 

Interval 

(cm-1) 

2000-1600 

1600-1100 

500-333 

333-250 

250-200 

200-0 

Spectral 

Interval 

(n) 

5 -6.3 

6. 3-9.1 

20-30 

30-40 

40-50 

50-=° 

5 (Si) 

gm 

2 . 2xl0 5 

6.6x104 

2.4xl0 5 

4.2x105 

1 . 1x10® 

a 0 (cm _1 ) 

0.10 

0.09 

0.09 

0.09 

0.09 

t (cm~l) 

1.0 

2.0 

2.0 

2.0 

2.0 

P'o 

0.6 

0.3 

0.3 

0.3 

0.3 

X(STP) + 

3.1xl0 3 

1.0x103 

3.8xl0 3 

6 . 6xl0 3 

1.7x104 


+ computed for 9 x 10“ 3 gm/cm 2 of HgO. 

2 1 1 Su 

+ Requirement for strong line approximation : X = — — — >1.6 
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Using fundamental data such as line widths, spacings and intensities, 
a check for the validity of the strong line approximation was 
carried out. This approximation was fitted to the data of Burch 
et al (Reference 3) for the 6.3d band and the data of Palmer 
(Reference 6) for the rotation band using water vapor amounts from 
10~3 to 10”2 gm/cm^. The resulting coefficients m and n used in 
the expression for the optical thickness in the spectral interval i, 
T i ~ (mj[Ut*) ni , are summarized in Table II, along with an updated 
list of CO 2 coefficients. The temperature dependence of the popu- 
lations of the rotational levels of the pure rotation band and of 
the vibrational level for the 6.3d band was neglected. The pressure 
and temperature dependence of the line half -widths are included. 

B. Vertical Distribution of Water Vapor 

The calculations described in Part I (Reference 1) incor- 
porate a constant mixing ratio for CO 2 . Since CO 2 condensation 
appears to be unlikely, this procedure is justified. However, the 
relative ease with which water vapor may condense for the pressures 
and temperatures considered, makes it desirable to represent the 
vertical distribution of its concentration in a variable fashion. 
This may be accomplished by specifying the concentration as a 
function of pressure, for example, as shown in Eq. (1) 


r q(P) , , H-H* . , P 

%— 1 - ( -TTF- > lQ g e 71 


The effects of a constant as well as a variable concen- 
tration will be discussed. 
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Table II - Absorption Coefficients 


i 

v (cm - -*-) 


in 

n 

Y 

1 

0-200 

2.6 

X 

10 3 

0.35 

0 

2 

200-250 

8.0 

X J 

jl02 

0.40 

0 

3 

250-335 

1.7 

X 

10 2 

0.38 

0 

4 

335-495 

1.2 

X 

10l 

0.35 

0 

5 

495-550 

7.8 

X 

10-5 

0.60 

0 

6 

550-625 

4.6 

X 

10“ 3 

0.42 

0 

7 

625-660 

9.0 

X 

io-i 

0.38 

0 

8 

660-720 

4.9 

X 

io-i 

0.42 

0 

9 

720-810 

9.5 

X 

io-4 

0.40 

0 

10 

810-880 

6.5 

X 

10" 5 

0.55 

2200 

11 

880-920 

3.2 

X 

io-7 

0.47 

2040 

12 

920-1000 

7.5 

X 

10- 6 

0.47 

2040 

13 

1000-1100 

1.6 

X 

10" 5 

0.51 

1950 

14 

1100-1600 

2.5 

X 

10 1 

0.44 

0 

15 

1600-2000 

1.0 

X 

iOl 

0.45 

0 

16 

2000-2600 

6.6 



0.50 

0 

17 

2600-8000 

1.1 

X 

10-3 

0.50 

0 

18 

1700-8000 (solar) 

4.9 

X 

10-7 

0.27 

0 


u* 

- co 2 

in 

cm atm 

u* 

- h 2 0 

in 

gm/cm 2 
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C. Convective Effects 


The results of Part I of this study showed in many cases 
a radiative gradient exceeding the adiabatic one, indicating a 
zone of instability. The radiative gradient is then not realized. 

The criterion for stability as formulated by K. Schwarzschild 
(See References 7 and 8, for example) is expressed in terms of the 
adiabatic gradient, and the gradient actually maintained. The 
atmosphere is stable against convection if 


( 2 ) 


(££) 

v dz' 


actual 




(dl) 

K dz J 


adiab 


This criterion and the flux continuity condition yield the inter- 
face between zones of radiative and convective transfer. For a 
homogeneous atmosphere these conditions are equivalent to continuity 
of the temperature and the temperature gradient at the desired 
interface . 


A modified form of Equation (2) was applied in the computer 
calculation using a method of successive approximations and taking 
the adiabatic profile as the initial temperature estimate. The 
atmosphere above a specific layer was permitted to adjust to radiative 
equilibrium to establish the gradient at this layer. Generally, 10 
iterations were found to be sufficient. At the same time, the adia- 
batic profile was assumed to exist between this layer and the cloud 
surface. All layers contribute to the radiative solution although 
the constant flux criterion is applied only to the layers above the 
convective zone. In most cases, the normalized gradient will be 
highly positive at the starting layer indicating that the desired 
boundary must occur at higher pressures. As the next step, the 
interface is adjusted toward the cloud surface and the radiative 
solution established above this new layer. This process continues 
until the boundary conditions are satisfied within a predetermined 
tolerance. After finding the boundary location, the solution in 
the radiative zone is iterated until the net flux is constant over 
the whole radiative zone. 

In addition to determining the tropopause height, the 
calculation permits an assessment of several other effects re- 
sulting from the convective transfer. At each point in the con- 
vective zone, the difference between the constant flux in the 
radiative regime and the radiative flux in the convective zone 
is the amount of energy transported by convection, rrF a . 
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In the atmosphere, this flux is produced by the vertical 
motion of mass elements, which may be written, neglecting geo- 
potential and kinetic energies, 

(3) ttF = pvc T 

a p 

In general, one half of the area element dA contains upward motion 
while the other half maintains downward motion, so that the net 
convective flux rrF a may then be expressed as 

(4) nF a = ^ p c p v(p-3*) 

The quantity is a characteristic length, v is the 
average velocity of the moving element, and (3-j3*) is the excess of 
the true gradient over the adiabatic gradient. 

To obtain a further relationship among the parameters the 
equation of motion of the atmospheric elements may be integrated. 
The equation of motion is expressed as the balance of the buoyant 
and the inertial forces acting on the element: 

(5) p* ^ = (p-p*) g, 

where p* is the density of moving element, and p the density of 
the environment. At any given time, the convective motions may 
be considered as the superposition of periodic, closed path loops 
carried out by the moving elements. To a first approximation, 

Eq. (5) may be integrated to yield an average velocity v, 

(6) v = ^ i rf } 1/2 • 

The average vertical velocity and the excess gradient may 
be reformulated by combining Eqs. (4) and (6) with the adiabatic 
gas law, to yield 

T-l g (h) ttF (h) 1/3 

(7) v(h) = — - and 

4T 2 P(h) 
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( 8 ) 


2/3 


B±|* =gH ( .1)V3 i_^a 

p p 2 

The quantities v and g* are expressed as functions of known quan- 
tities with the exception of the convective flux, rrF a , and the 
mixing length, 

The convective flux was obtained directly from the computer 
solution. The mixing length, , is more difficult to obtain. It 
may be defined as the characteristic length over which the linear 
momentum, pv, is conserved, approximately. In a vertical column of 
gas the density p and the velocity v of the rising elements steadily 
diminish with height. Hence, an outflow of material must occur at 
the sides of the column. The distance over which 


(9) 


P2 V 2 ^ \ Pl v l 


is defined here as the mixing length and may be expressed more 
conveniently as 


( 10 ) 


1 _ 2 d(pv) 

pv dh 


By an iterative process involving Equations (7) and (10) and the 
computed convective flux, the coefficients for a Taylor expansion 
of the mixing length are obtained, 

Ui) (h) = \ *<h) - g H * 2 < h > > 

where x is the distance from the tropopause and H the density 
scale height of an adiabatic atmosphere. The distance x is assumed 
to be small compared to the density scale height. The mixing length 
is slightly less than one half of the distance from a particular 
layer to the tropopause. 

With this result, the superadiabatic gradient and the 
average vertical velocity may be computed. The gradient was found 
to exceed the adiabatic gradient by one degree per km at most, so 
that the use of the adiabatic profile in the computation appears 
justified. The desired vertical velocities were then computed 
using the mixing length according to Eq. (11) . 
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One may then proceed a step further and determine the 
maximum particle size which could be supported by these velocities 
on the assumption of Stokes’ law: 

: 1/2 

(12) d = (i|3L£) 

p« s 


For the rough estimates of particle sizes presented later, more 
sophisticated forms of Stokes’ law such as used in Ryan’s study on 
Mars (Reference 9) are probably not warranted. 

D. Temperature Calculation 

The solution for the temperature distribution described 
previously (Reference l) requires the specific temperature which 
balances the emission of each layer with its absorption. A relation 
ship between the emission of a layer and its temperature was estab- 
lished by computing an array of values of emission, optical 
thickness, and temperature. When the absorption was computed, the 
value of emission identical to the absorption with the appropriate 
value of optical thickness was used to inspect the previously 
computed table for the value of temperature which satisfied the 
above conditions. This procedure sometimes yielded slightly 
different values of optical thickness in each side of the energy 
balance equation for those cases which showed very steep gradients. 

A new approach was therefore adopted in which the value 
of optical thickness for a layer, based on the previously iterated 
temperature distribution is used in the emission as well as the 
absorption term, thus eliminating the small discrepancy. 

III. EFFECTS OF INDIVIDUAL PARAMETERS 


The maximum number of parameters used in any analysis is 
determined by the ability to represent the physical situation and 
by the number of independent observational data which are available 
for interpretation. For this study of the Venus atmosphere, the 
following parameters were selected. 
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Cloud Top Pressure 

- p c 

Cloud Top Temperature 

- T c 

Cloud Top Spectral Emissivity 

e ic 

Cloud Top Reflectivity 
(near infrared, l-5|~i) 

- r c 

Carbon Dioxide Concentration 

- <ico 2 

Water Vapor Concentration 

- <ih 2 o 

Solar Zenith Angle 

- c 


The effects of a variation of specified combinations of 
these parameters may be studied on temperature distributions, thermal 
emission spectra, and limb functions. In the following discussions, 
the effects of the above parameters on each of these results will be 
presented separately. 

A. Effects on Temperature Distribution 

The results of a previous parametric analysis (Reference 
1) applied to the equilibrium temperature distributions are summarized 
in Figures 1 and 2. In Figure 1, the strong effect of the near- 
infrared solar heating is evident from the change of temperatures 
with the solar zenith angle. The large temperature differences 
found between night and day are understood by a larger solar constant, 
a relatively high CO 2 concentration, and a high cloud top reflectivity 

Figure 2 illustrates the influence of the CO 2 concentration 
on the equilibrium temperatures for two cloud top pressures. The CO 2 
effects differ somewhat for night and day. On the night side tem- 
peratures depend little on the CO 2 concentration, except in the upper- 
most layers, which are able to cool more effectively for the higher 
CO 2 concentrations. The illuminated hemisphere shows larger effects, 
especially near the cloud surface where the opacity for solar 
radiation is greatest, thus increasing the solar heating in these 
layers for higher CO 2 amounts. The strong influence of the cloud 
top pressure is also evident. 

The influence of the cloud top emissivity op the tem- 
perature distributions is not illustrated. An emissivity of 0.8. 
showed negligible influence on the shape of the equilibrium tem- 
perature distributions. Similarly, a 10% variation in the adopted 
value of 0.6 for the near infrared cloud top reflectivity shows a 
minor effect and is not shown. 

Water vapor influences the equilibrium temperatures 
distinctly. This is shown in Figure 3 for a cloud top pressure 
of 0.3 atm, for noon and night conditions. Several concen- 
trations and their respective saturation curves (Reference 10) 
are illustrated. 
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Also shown for the same total amounts of H 2 O under noon 
time conditions, are two additional temperature distributions 
using water vapor scale heights of 1/3 and 2/3 the pressure scale 
height. Such variations allow an effective balance of H 2 O cooling 
and CO 2 heating, yielding a nearly isothermal atmosphere. At the 
same time, the temperature distribution in the upper regions 
depart further away from the saturation curves. Tropopause heights 
are also indicated, when they appear. 

The effects of water vapor on the equilibrium temperature 
distributions may be summarized as follows. 

1. Large amounts of water vapor introduce sub- 
stantial cooling of the sunlight hemisphere. 

This is understood from the atmospheric 
opacity in the rotational band, where approxi- 
mately one-half of the energy spectrum lies.. 

2. A growth in the convective zone results with 
the addition of water. Again, the cooling 
influence of water vapor increases the radiative 
gradient and the convective flux, thus 
extending the zone to greater altitudes. 

3. The condensation properties of water vapor 
coupled with the convective instability imply 
that the night side cloud tops would grow to 
higher altitudes than the day time clouds, 
other parameters being held constant. On 
the illuminated hemisphere condensation may 
or may not take place depending on the water 
vapor concentration and its vertical dis- 
tribution. In either hemisphere, the latent 
heat released or absorbed in the vapor-ice 
phase transitions would yield a substantial 
influence on the resulting cloud top temperature 
as well as the vertical distribution adjacent 

to the cloud tops. 

As discussed previously, an assessment of the effect 
of the atmospheric parameter on the relevant quantities of the 
convective motion is possible. The dependence of the tropopause 
height, the magnitude and the vertical distribution of the con- 
vective flux, and the associated vertical gas velocities may be 
studied. Figure 4 shows each of these items for one particular 
combination of atmospheric parameters. 

The situation illustrated is for an 0.3 atm cloud top 
pressure level under night time conditions and unit emissivity. 

The convective flux, the mixing length, the mean vertical velocity, 
and maximum particle size are shown as a function of the distance 
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from the tropopause, In all cases studied the convective flux 
increased with the square of the distance from the tropopause. 

Since the analysis is based upon a mixing length theory and other 
simplifications, the results for the mean velocities and particle 
sizes are only approximate. 

A study of several sets of computations shows that the 
magnitude of the convective flux increases as the cloud top 
pressure is increased, primarily as a result of the increased 
atmospheric density. The corresponding mean vertical velocities 
and maximum particle diameter decrease in magnitude, however. 

C. Effects on Thermal Emission Spectra 

The effects of the atmospheric parameters on the thermal 
emission spectra will now be discussed. The most influential 
parameter is the solar radiation. For the case illustrated in 
Figure 5 where the brightness temperature is plotted as a function 
of wavenumber, with the wavenumber scale arranged in equal energy 
increments for a 225°K blackbody, solar heating reverts the 15 m. CO 2 
band from pure absorption to a partially absorption and partially 
emission spectrum. The most dominant feature is the 15 m CO 2 band, 
although the rotation band of H 2 O appears moderately strong. The 
8-13U region appears rather uniform, except for subsolar conditions 
where the temperature dependence of the CO 2 bands, combined with a 
large temperature inversion produce moderate features. 

The effect of cloud top pressure is also pronounced 
(Figure 6). For low pressures, very little structure appears in 
the spectrum, due to the small total opacity of the atmosphere. 

At higher pressures spectral features stand out strongly, even for 
small CO 2 concentrations. Figure 7 shows the effects of varying 
the CO 2 concentration under night time conditions. As expected, 
the spectral features become more intense as the opacity of the 
atmosphere increases. 

The effects of cloud top emissivity are shown in Figure 8. 
For this study, the emissivity beyond 13M was taken as unity but 
was permitted to assume some specified value at shorter wave- 
lengths. The combination of surface temperature and emissivity 
was adjusted to maintain about 225°K, close to the observed value 
of the brightness temperature in the atmospheric window. The 
effects in the 15 m band and the rotational band are moderate. A 
small increase in the surface temperature with a constant emissi- 
vity elevates the entire spectrum, keeping the shape of the spectrum 
essentially intact. 


12 



Figure 9 illustrates the effects of varying the H 2 O 
concentration. The primary influence appears in the rotational 
band for large concentrations. 

C. Effects on the Limb Functions 

The computed limb functions allow further studies of the 
model atmospheres. The limb functions depend primarily on the 
optical thickness of the atmosphere. In the absence of strong 
scattering effects, the variations in the limb function depend on 
the temperature variations over the surfaces of constant optical 
thickness . 


Two of the eighteen spectral intervals are illustrated in 
Figure 10. These are the 1064 cm~l band and the 625-640 cm -1 
region and they represent weakly and strongly absorbing intervals, 
respectively. The limb function variation is determined by the 
magnitude of the temperature inversion, and therefore by solar 
heating, and the opacity of the spectral interval. The effect of 
solar heating on the weak interval is well illustrated while the 
COg influence is shown best by the strongly absorbing interval. 

Figure 11 shows the effect of cloud top pressure on the 
limb functions, the variations being larger for the weakly (495-550 
cm”l) absorbing interval. The H 2 O variation is small, although a 
slight reversal of the pressure effect is seen in the strongly 
absorbing interval (625-660 cm-1) in Figure 11. The effect of 
varying the cloud top emissivity is negligible for the intervals 
discussed and is not illustrated. 

IV. APPLICATION OF RESULTS TO A COMPARISON WITH SOME VENUS 

OBSEWATTcKs 

Several observations are available which may provide a test 
of these calculations. These are considered in the following 
order: 

A. Various observations relating to the presence of 
water vapor. (References 11, 12, 13, 14, 15, 16, 17) 

B. Limb function observations in the 8-13n interval by 
Sinton and Strong (Reference 18) , by Sinton 
(Reference 19) and the more recent thermal maps 
obtained by Murray, Wildey, and Westphal (Reference 
20 ) . 

C. The thermal emission spectra measured in the 8-13M 
interval by Sinton and Strong (Reference 18) . 
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A. Water Vapor Measurements 

A knowledge of the presence of water vapor in the Venus 
atmosphere is very significant. In moderate to large amounts, water 
influences the atmospheric structure from a radiative as well as a 
dynamical point of view. The question of cloud composition is 
also associated with the presence of water in the vapor state. 
Finally, the presence of water is an important indication of the 
historical development of the planet and its atmosphere. The dis- 
cussion presented here will only be concerned with the influence of 
water vapor on the atmospheric structure and its relation to 
cloud phenomena. 

The existence of water vapor in the Venus atmosphere was 
first suggested by Strong et al’s 1959 balloon measurements 
(Reference 21) . These were repeated in 1964 with improved apparatus 
and again yielded positive results (Reference 13) . Similar measure- 
ments by Dollfus (Reference 17) using the 1.38|j region were also 
positive. In both cases the actual quantities of water vapor 
deduced require considerable interpretation since the values depend 
heavily on the pressure level of the reflecting layer. On the other 
hand, Spinrad's (Reference 11) spectroscopic analysis of the .87p. 
bands has produced 10 _ 6 f or the mass mixing ratio, a value leading 
to amounts below Bottema et al's and Dollfus’s amounts, unless the 
reflecting cloud layer is placed at relatively high pressures 
(>1 atm) . 


Another observation indicating the presence of water vapor 
is the more recent balloon measurements of Bottema, Plummer, Strong, 
and Zander (Reference 14) which have shown that the reflectivity of 
the cloud layer in the Venus atmosphere to be similar to that of 
ice crystals. Nevertheless, high resolution spectra obtained by 
Kuiper (Reference 15) and by Moroz (Reference 16) near 2|i do not 
reveal the characteristic ice crystal absorption. No attempt will 
be made at reconciling these results; instead an analysis pre- 
senting both viewpoints will be considered and illustrated here. 

Since the deduced amounts of water vapor depend strongly 
on the cloud top pressure level, the limits usually quoted for this 
level in the literature (Reference 22, for example) (0.09 -0.60;atm) 
will be illustrated. Table III summarizes the water vapor con- 
centrations deduced for each pressure level, for two of the investi- 
gations cited. The results of the calculations using the values 
listed in Table III are illustrated in Figures 12, 13 and 14. 
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TABLE III - HoO 


Reflecting Layer Pressure Level 

Spinrad(1963) Volume Concentration 

Bottema et al (1965) Volume 
Concentration 

Saturation Volume Concentration 
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Consider first a cloud layer consisting of ice crystals.' 
Under these conditions, the observed amount of water must provide 
saturation at the cloud top. For 0.09 atm (Figure 12) and a cloud 
top emissivity of unity, the large degree of cooling generated by 
the emission in the H 2 O rotational band offsets the solar heating 
due to near infrared absorption by CO 2 . Comparison with the con- 
densation curve implies a super-saturated atmosphere even under 
noon-time conditions and hence is incompatible with this pressure 
level. On the night side, the condensation coupled with convective 
motions would elevate the cloud top to higher altitudes. 

In Figure 13, the reflecting cloud layer is placed at 0.60 
atm and the results show that the solar heating produced by the near 
infrared absorption of CO 2 dominates the H 2 O cooling. The tem- 
perature distributions do not allow saturation at this pressure 
level. Hence, the cloud layey, if composed of ice, must be placed 
at a higher elevation, until saturation is possible at the cloud 
tops. This occurs at approximately 0.30 atm, although a more 
definite level is not offered since the absorption law (1.13ii band) 
is intermediate between weak and strong for this pressure range 
and water amount. This interpretation implies that the observed 
water pertains essentially to the clear atmosphere above the clouds 
and not to deeper layers. 

The above analysis may represent an over-simplified inter- 
pretation since other atmospheric conditions are capable of 
providing alternative points of view which might still be consistent 
with the 0.09 atm level. Among these possible conditions are: 

1. The cloud top emissivity may be less than unity and 
therefore the true temperature in the vicinity of the cloud tops 
may be higher than the observed brightness temperature of about 225°K. 
Cloud emissivities near unity are characteristic of thick water or 
ice crystal clouds. Lower values of emissivity correspond to either 
nonaqueous or semi-transparent cirrus-type clouds. The reduced 
apparent brightness temperature which results in such cases corres- 
ponds to higher values of true cloud top temperature. The corres- 
ponding equilibrium temperature distributions may then be shifted 
beyond the saturation curves and condensation thus avoided. 

The calculations using a cloud emissivity of 0.5 with 
a cloud top temperature of 255°K placed at 0.09 atm is shown in 
Figure 14. Condensation would again occur, and since 0.5 is 
probably a lower realistic limit for the cloud emissivity, this 
pressure level may be rejected. 


16 



2. The H 2 Q 1 measurements may have actually "penetrated" 
through the clouds to deeper layers, so that the water vapor 
concentrations shown in Table III should be considered upper 
limits. If water amounts which could avoid condensation are desired, 
then the observation would be required to penetrate one or two 
scale heights into or below the clouds. Although this requirement 

is rather severe for a single reflecting cloud layer, it is 
difficult to reject categorically. 

3. The water vapor may exist in a nonuniform vertical 
distribution in which case the total amount of H 2 O deduced would 

be modified. Also shown in Figures 12 and 14 are temperatures com- 
puted with a very nonuniform vertical H 2 O distribution, keeping the 
total amount the same as in the uniform distribution case. For a 
cloud emissivity of unity, the previous conclusions are not signi- 
ficantly altered. If the cloud emissivity is about 0.5, then con- 
densation for noon time conditions could be partially avoided, 
although larger solar zenith angles would still exhibit condensation. 
Again, this constraint on water clouds is severe, and the vertical 
distribution used rather extreme. Condensation appears to be diffi- 
cult to avoid. Hence, the third alternative may be discarded. 

The net conclusion of this part of the analysis, with some reser- 
vation, is that the reflecting ice crystal cloud layer should exist 
at the 0.20-0.40 atm pressure level with an assumed cloud top 
temperature of 225°K. 

A similar analysis may now be considered for nonaqueous 
clouds. In this case, the atmosphere must everywhere lie well 
below saturation. A typical value that may be considered for the 
mass mixing ratio of water vapor is Spinrad’s value of 10~6. For an 
emissivity of unity, condensation is not possible under daytime 
conditions for either the 0.09 or 0.60 atm cases as shown in 
Figures 12 and 13. For both pressure levels, condensation and 
cloud formation would be likely on the dark hemisphere although 
an emissivity as low as 0.5 may restrict the condensation on the night 
side to thin haze layers in the upper regions. Concentrations less 
than 10“8 or emissivities of about one-half would be necessary to 
avoid H 2 O clouds entirely on the night side. Either requirement 
would be consistent with nonaqueous clouds for most pressure levels. 
However, two limits on the cloud top pressure may be derived on the 
assumption of nonaqueous clouds. If a total water content of 
9 x 10 “3 gm/cm^ (Bottema et al) (Reference 13) is assumed to exist 
above the clouds, with a cloud top temperature of 225°K, then cloud 
top pressures in excess of 0.4 atm are required to avoid saturation 
at the cloud tops. An upper limit may be determined with 
Spinrad’s mixing ratio of 10“”. An examination of Figure 13 
reveals that a temperature of 225°K is not achieved with the 
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saturation curve for pressures less than 1 atmosphere. In fact, a 
pressure of 31 atmospheres is necessary to reach a temperature of 
225°K. If lower emissivities and higher temperatures are 
appropriate, then this pressure range would be shifted to higher 
values. Further conclusions on the cloud level and other atmos- 
pheric parameters will be considered in the subsequent sections 
dealing with other observations. 

B. Limb Function Observations 

The detailed thermal maps produced by Murray, Wildey and 
Westphal (Reference 20) and the limb functions obtained by Sinton 
and Strong (Reference 18) and by Sinton (Reference 19) provide 
another set of observational data to compare with the computed 
results. The latter measurements indicate a darkening proportional 
to the square root of the view angle cosine, while Murray et al's data 
show a somewhat weaker dependence. 

The results are compared in Figure 15 in which the par- 
tially illuminated hemisphere is observed to yield a darkening not 
generally different from the night side of the planet. However, 
the calculations using intensities averaged over the six spectral 
intervals comprising this portion of the spectrum indicate a sub- 
stantial brightening in the illuminated hemisphere. The illuminated 
hemisphere is represented by a composite of two zones, whose average 
solar zenith angles are 60° and 75°. These two zones used in con- 
junction with a dark side computation are sufficient to represent 
the planetary phase during the observations of Murray et al . Three 
situations are compared with the observed result. The light solid 
curve results from a calculation placing the cloud tops at 1.0 atm. 
The dashed and dotted curves result from an 0.3 atm pressure level 
under various assumptions. The dashed case results from differences 
in the computed temperature distributions over the disk and show 
little variation with some asymmetry in the limb function, while the 
1.0 atm case shows appreciable night side darkening with consider- 
able asymmetry. The dotted case is discussed later. A darkening 
comparable to the observed amount can be obtained with a more 
opaque atmosphere (P c > 1 atm). Such an opaque atmosphere, how- 
ever, yields an even greater degree of brightening on the illumi- 
nated hemisphere, which is enhanced and the asymmetry increased in 
the computed limb functions as the planetary phase is increased. 
Sinton ’s observations clearly illustrate that little change, if 
any, occurs in the shape of the limb function under various phase 
angles. Hence, the solar heating and its related brightening 
must be minimized or its effects eliminated to yield accordance 
with the observations since an explanation of the observed darkening 
cannot be made on the basis of an absorbing atmosphere alone. The 
disparity between the observations and calculations suggest the 
following considerations as means for reconcilation . 
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1. A very low CC >2 concentration (qcOo <<: ^.05) may be 
necessary to minimize the effects of solar heating and its subse- 
quent contribution to the limb brightening. However, such low 
concentrations, whatever their value, must still be capable of 
yielding the observed equivalent widths of the CO 2 bands in the 
red and near infrared portions of the spectrum. Moreover, the 
atmosphere becomes so transparent, unless large values of cloud 
top pressures are used, that the requisite absorption is not 
available to yield the darkening function. 

2. The clouds may be responsible for the observed limb 
darkening. Scattering by cloud particles, a large negative tem- 
perature gradient within the cloud, the vertical cloud structure, 
or some combination of these items may contribute to the observed 
limb darkening. A number of authors have considered these 
possibilities (References 23, 24,25,26 and 27 for example) and have 
demonstrated that each is capable, alone or with some realistic 
combination, of explaining the observed darkening. 

3. If larger CO 2 amounts do exist, then the computed 
brightening in the sunlit hemisphere may be removed by dynamical 
processes. Global convection currents giving rise to adiabatic 
cooling in the sunlit hemisphere and an adiabatic heating on the 
night side could yield approximately the same temperature dis- 
tribution above the clouds for day and night conditions. Solar 
heating may be sufficient to act as a driving force for a general 
atmospheric circulation. The maximum effects permitted by such 

a circulation, with the observed amounts of H 2 O, would yield a 
temperature distribution somewhat above the corresponding 
saturation curve for H 2 O. Although this distribution is rather 
arbitrary, it is a convenient and physically plausible one to use. 
The net result of such circulation may be simulated with these 
computations and the result is shown as the dotted curve in 
Figure 15. The degree of darkening still falls short of the 
required amounts. High pressures, coupled with such circulation 
would be required, but saturation cannot be achieved at these 
levels with the observed amounts of H 2 O. 

The conclusions from this part of the analysis hinge on 
the presence of several possible physical conditions or phenomena. 
If the clouds consist of ice crystals then somewhat lower CO 2 
concentrations, scattering by the cloud particles, or atmospheric 
motions to remove the solar heating effects, or a proper com- 
bination of all three items are required to yield the observed 
darkening. 
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If the clouds are not water vapor, then even lower CO 2 
concentrations would be necessary to reduce the solar heating, but 
then the observed darkening becomes more difficult to achieve on 
the night side. The net result remains the same, in either case, 
that the clouds must be responsible for producing the observed 
limb darkening. 

C. The rmal Emission Spectrum of Sin ton and Strong 

Measurements of the emission spectrum of the planet Venus 
were reported by Sinton and Strong in 1960 (Reference 18). The 
total disk of the planet was used to record the spectrum from 8-13|i. 
In order to simulate this observation with these calculations, the 
planet was divided into several zones of varying solar zenith angle. 
In addition, identical zones of aspect or view angle, rotated 
ninety degrees from the zenith angle zones toward the observer, are 
superposed on the planet. The resulting grid is shown in Figure 16, 
where the mean solar zenith angle and mean aspect angles for each 
zone are indicated. Symmetry permits consideration of only one 
half the observed hemisphere. The computed spectra were synthesized 
by taking the specific intensity from each of the six areas and 
weighting the values by the relative projected area of each zone. 
These are then summed to yield a flux value emitted by the illumi- 
nated hemisphere. This is repeated for each spectral interval over 
the range 8-13ia, using intervals 40-100 cm"-*- wide. These fluxes 
are then averaged with their corresponding spectral fluxes from 
the night hemisphere to yield the final value. A refinement of 
the procedure using greater spatial resolution yielded no signi- 
ficant differences. 

Ideally, a comparison yielding agreement of the computed 
spectra with the observed one, is desirable over the whole 8-13ii 
interval. At best however, only a qualitative agreement could be 
achieved in two ways. The overall shape of the spectrum, neglecting 
the 11.2ia feature can be reasonably computed. Alternatively, a 
comparison can also be attained on the basis of the 11.2|j feature 
alone. Several combinations of atmospheric parameters were con- 
sidered and synthesized, and are shown in Figures 17 and 18, where 
the brightness temperature is plotted as a function of wavelength. 
The effect of cloud top pressure is shown in Figure 17, assuming a 
cloud top emissivity of unity. The overall shape of the spectrum 
is fairly well reproduced, although as may be expected, no feature 
at 11.2ia is apparent. A slight increase of cloud top temperature 
will increase the general level of the spectrum to yield a better 
agreement, without altering its shape. It may be noted that an 
increase in the pressure level yields colder temperatures in the 
12-13U region. All three pressure levels appear to yield fair 
agreement . 
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Figure 18 illustrates an attempt at duplicating the 11.2 m 
feature, without the presence of an additional atmospheric absorber 
or a peculiar cloud top spectral emissivity. A cloud top emissi- 
vity of one-half is shown for several pressure levels. Good agree- 
ment in the 10-13M region is obtained with a cloud top pressure of 
0.3-0. 5 atm. Again, the general level may be increased to improve 
the match by a slight increase in cloud top temperature without 
destroying the shape. Emissivities smaller than 0.5 yield excess- 
ively high temperatures in the 8-10|i region while values larger than 
0.7 are insufficient to yield the 11 m feature, which originates from 
a combination o,f the transparency of this spectral interval and the 
decreased emissivity. The 0.3-0. 5 atm pressure level appears to 
yield the best fit for this portion of the spectrum. 

At this point it may be mentioned that a more detailed or 
extensive study is not justified until more observational studies 
confirm the spectrum and its 11.2 m feature, preferably under vary- 
ing spectral resolution, if possible. Pending such confirmation, 
however, it is felt that the general trends introduced into the 
spectra by various factors are still of interest and may contribute 
to the number of qualitative deductions permitted by such an 
analysis . 


A factor necessary to consider is the influence of the 
observed limb darkening on the spectra. This effect may be 
considered by applying the observed limb function to the synthesized 
spectra, such that the contribution from the outer zones of aspect 
angle is reduced. In effect, the weighting factors for the outer 
zones were decreased according to the observed limb darkening. Dis- 
cernable differences of spectra computed with and without this 
alteration were not found. 

A second factor which influences the spectra is the pre- 
sence or absence of solar heating on the sunlit side of the planet 
arising from near infrared CO 2 absorption in the atmosphere. If 
radiative equilibrium alone prevails then the sunlit side of the 
planet should be warmer than the night side and the resulting 
spectra are those shown in Figures 17 and 18. If, however, planetary 
convection predominates over the radiative equilibrium, nearly 
uniform temperatures may be expected on the day and night sides of 
the planet and result in correspondingly different spectra. This 
result may be simulated in the computed spectra by examining the 
spectra resulting from the use of a single temperature distri- 
bution over the planet. An arbitrary but physically reasonable 
temperature distribution to use under such assumptions in the pre- 
sence of water vapor, would be the saturation curve corresponding 
to the observed amount of water. An examination of the spectra 
resulting from the use of the H 2 O saturation curves as the tem- 
perature distribution in both hemispheres yields a flat blackbody 
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spectrum of about 225°K in the 8-13p. region. The use of a steeper 
gradient gives nearly identical results, with only a slight drop 
in the 12-13P temperatures. Hence, a composite spectrum as shown 
in Figures 17 and 18, which includes the differences in the day 
and night hemispheres, appears to be necessary to yield a close 
agreement in the overall shape of the spectrum. 

The effect of the accuracy of the absorption coefficients 
may also be discussed presently. A study showed that a 20% 
change in the coefficients yielded no significant difference in 
the shape of the computed spectra. 

Further conclusions based on a spectral comparison alone, 
require a discussion of the cloud composition. Cloud emissivities 
near unity, which are characteristic of opaque ice crystal clouds, 
are capable of yielding computed spectra in fair agreement with 
the general shape of the observed spectrum, using a cloud top 
pressure in the range 0.1-0. 3 atm. The 11.2|u feature is not 
apparent under these conditions. On the other hand, cloud emissi- 
vities in the range 0.5-0. 7 correspond to nonaqueous clouds and 
provide computed spectra yielding the llu feature, but with some- 
what poorer agreement in the 8-9u interval, using cloud top 
pressures from 0.3-0. 5 atm. Spectral emissivities peculiar to 
other cloud substances, in combination with the other atmospheric 
parameters cannot be ruled out, and may yield better agreement. 

An interesting possibility is a combination of ice crystallized 
onto particles of CaCOg or Na 2 CC> 3 „ Either of the latter two 
substances show strong features at 11. 2u and may be possible 
cloud particle candidates (Reference 28, 29) . 

These conclusions are at best qualitative and suggestive 
only. A factor not considered here which may play as great a role 
in determining the observed spectrum as the composition, is the 
vertical ‘cloud structure. Further analysis must await repeated 
observation of the spectrum, preferably carried out over various 
portions of the disk, isolating the sunlit and dark hemispheres, 
if possible. 

V. SUMMARY AND CONCLUSIONS 

An attempt has been made to illustrate the important effects of 
several parameters on the atmospheric structure in the Venus atmos- 
phere above the cloud layer. The results indicate that the cloud 
top pressure, water content and solar heating by CO 2 absorption 
have a strong influence on the atmospheric structure. A comparison 
of these calculations with several observational data provides a 
severe test. 
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Within the limitations of the calculations and the assumption 
of a 225°K brightness temperature for the cloud tops and a 5% CO 2 
concentration, the results of such a comparison lead to the 
following conclusions: 

1. A study of the water vapor abundance yields a cloud 
top pressure of ~0.3 atm for ice crystal clouds, 
and pressures greater than 0.4 atm for nonaqueous 
clouds . 

2. Studies of the thermal emission spectrum show good 
agreement with the observations, over the pressure 
range, 0. 1-0.8 atm, for unit cloud emissivities and 
0.3-0. 5 atm for smaller emissivities (~0.5) . 

3. The limitations on the cloud top pressure level may 
be further restricted by merging the above results. 

A pressure range of 0,2-0. 4 atm appears most 
probable for ice crystal clouds while values from 
0.4-0. 5 atm appear reasonable for nonaqueous clouds. 

These limitations on the cloud top pressure are based 
on a constant spectral emissivity. However, the clouds 
themselves may be totally responsible for the observed 
spectrum, in which case the use of an appropriate, 
variable spectral emissivity would probably yield 

some adjustment in these limitations. 

4. An analysis of the limb function calculations yields 
no further limitations on the cloud top pressure but 
yields the conclusion that the clouds are responsible 
for the observed limb darkening. Earlier studies, 
not considered here, of the planetary heat balance 
yielded similar inconclusive results for the cloud 

top pressure, within the limits of the observed albedo. 

In order to improve the precision of the results and to clarify 
the importance of various atmospheric phenomena and other factors, 
the following observations are suggested. 

a. A continuation of the detailed thermal mapping 
as a function of planetary phase. 

b. A repetition of the 8-13|J. spectrum using distinctly 
sunlit or night portions instead of the total planet. 

c. If possible, measurements of the total thermal 
spectrum including the 15p. C0 2 band, again re- 
solving day or night portions of the planet. 

d. Ground-based measurements of the near IR C0 2 bands 
would also be desirable as a function of solar zenith 
angle over the planet. 
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e. High spatial and spectral resolution studies of 
the photographic IR spectrum of Venus. Possible 
variations with time and position over the sunlit 
disk of the planet, as a function of phase, are 
needed to provide a better physical picture of 
the cloud structure, as well as a clearer 
indication of the relative roles of scattering 
and absorption. High resolution spectra with the 
slijt placed over various zones of equal solar 
zenith angle may provide some indication of 
variations in the rotational temperature over the 
sunlit part of the disk. Similar spectra with 
the slit adjusted perpendicular to the terminator 
may yield some data on variations in the cloud 
top elevation with solar zenith angle, if present. 

As may be seen from the above conclusions, the incomplete 
agreement of the calculations with certain observations suggests 
that the simple reflecting cloud layer model may be inadequate. 
Hence, further theoretical studies in developing more realistic 
cloud models and their subsequent influence on the atmospheric 
structure are needed. 
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Figure 1 —Radiative equilibrium temperatures for several solar zenith angles. Comparison with the 
adiabatic profile indicates instability for large zenith angles. 







Figure 2-Radiotive equilibrium temperatures for both night and day, and for various pressure levels 
and C02 concentrations. Condensation of CO2 is unlikely for the conditions shown. The adiabatic 
profile (5% CO2) is shown. 



Figure 3 — Radiative and convective equilibrium temperatures for several water vapor concentrations and 
water vapor scale heights. Also indicated are water vapor saturation curves and tropopause heights. 





Figure 4— Convective flux, mean vertical gas velocity, mixing length, density and 
maximum diameter of particles which can be suspended. 
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Figure 5— Thermal emission spectra for 
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■Thermal emission spectra for several cloud top 
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Figure 7 — Thermal emission spectra for several CO 2 concentrations in the dark hemisphere. 



T(°K) €=0.5 1 € c = 0.5, NOON 



36 


Figure 8— Thermal emission spectra for two cloud emissivities for the dark and the illuminated hemi 
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-Thermal emission spectra for several water vapor concentrations in the dark hemisphi 
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Figure 10— Planetary limb functions in a weakly and a strongly absorbing spectral interval. 
Various CO 2 concentrations are shown for both hemispheres. 
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Figure 1 2— Equilibrium temperatures for an .09 atm cloud top pressure, emissivity of one, and water 
vapor concentrations corresponding to the data of Bottema, Plummer and Strong as well as to the 
data of Spinrad. 
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measurements of planetary spectra. 
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observed spectrum of Sinton and Strong. 
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Figure 18— Synthesized spectra in the 8-13^ window for several cloud top pressures 



